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by scattering second sound off quantized vortices in superfluid Helium is sug- 
gested. Theoretical calculations of the relative amplitude of the scattered second 
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presented. The relevant estimates show that an experimental verification of the 
method is feasible. Moreover, it can even be used for the detection of a single 
quantum vortex. 
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Quantum vorticity in a superfluid, its properties, nucleation and dynamics, 
has been an important and long standing subject in a quantum fluid physics 
for several decades 0. The main experimental tool from the very beginning to 
investigate quantum vorticity has been second sound0. Ever since the first ex- 
periments by Vinen0] the only method which has been used in this respect, is 
the measurement of the variation and anisotropy of second sound attenuation due 
to the presence of quantum vorticity Q. This has turned out to be a very useful 
and powerful tool. However, the fact that second sound is a radiative hydrody- 
namic mode suggests alternative methods to probe quantum vorticity, by analogy 
with the methods that have been developed in classical hydrodynamics in recent 
years. One of them is to measure the phase difference of a dislocated, or splitted, 
wave front, a phenomenon that is analogous to the Aharonov-Bohm effect This 
method has been used to probe a bathtub vortex with water surface waves0. It 
was also used recently to measure the quantized vorticity of a single vortex in 
superfluid helium]?]]. Another approach, which is the topic of this Letter, is to 
measure scattering of sound waves by vorticity ||. This method was successfully 
used recently to measure vorticity in the vortex array of the Karman street gener- 
ated by flow past a cylinder at low Reynolds number 0. In this Letter we present 
calculations and numerical estimates for the scattering of second sound waves due 
to quantum vorticity. Calculations are valid for any number of vortices both in 
two and three dimensions, and in three they may be of arbitrary shape. The es- 
timates convincingly show that the method discussed is very sensitive and allows 
for the detection of a single quantum vortex. This method has the advantages of 
being practically unperturbative and fast, so it can be used as a powerful probe 
to study nucleation of quantum vortices or vortex rings in superfluid helium. 

A single, straight, isolated quantum vortex consists of a core of radius a, out- 
side of which the flow is potential, with the velocity of the superfluid component 
given by 

v s = V0 

where = T6/2tt, with 9 the polar angle around the point vortex, and V = h/m 
the quantum of circulation. It is important to note that this velocity potential 
is multivalued. At this level the normal flow is insignificant. Next, a second 
sound wave is sent onto this vortex. The wave motion will induce an oscillatory 
motion of the vortex core which will in turn radiate second sound. This is the 
basic mechanism for the scattering of second sound by a quantum vortex. The 
scattering of first sound by a quantum vortex in two dimensions was calculated by 
Pitaevski |10|| in a Born approximation, and by Fetter |nj in an exact form. Both 
used a hydrodynamic approximation. This issue has been recently readdressed 
within the framework of a many body wave function descrition of elementary 
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exitations in UT2 . 



In order to obtain a first approximation to this problem from a formal point 
of view, it is best to fix one's attention upon a potential velocity field, soon to be 
identified with the superfluid velocity, v s = V0, whose scalar potential </> obeys 

° 2 ° clV 2 4> = (1) 



dt 2 



with C2 the speed of (second) sound waves, for time scales such that viscous 
and thermal conduction effects can be neglected. Vortices are accounted for by 
multivalued velocity potentials 0: in two dimensions, a polar angle centered at Xq 
gives a point vortex in that location; in three dimensions the solid angle subtended 
at a given point by a closed curve gives a vortex loop along that curve. Nonlinear 
terms can also be neglected when the second sound frequency is large compared 
to the time scale for vortex flow on a length scale comparable to wavelength A. 
Also if A ^> a, the vortex can be considered as a point in two dimensions, or as a 
filament in three. The problem of finding the second sound generated by a vortex 
filament undergoing prescribed motion is thus reduced to finding the appropriate 
multivalued solutions of Eqn. (|l|). The necessary theoretical framework to handle 



such problems has been developed by one of usp3[ in the context of classical 
hydrodynamics, and the relevant formulae can be simply translated. The main 
result to be used is that the radial component of the far field velocity, v r s1 generated 
by a vortex filament parametrized by X(a, t), where a is a Lagrangean parameter 
labeling points along the filament: 

where an overdot means time derivative, a prime means derivative with respect 
to a and the right- hand- side is supposed to be evaluated at the retarded time 
t ret = t — \r\/c2- If a plane wave of amplitude v^, wave vector k Q and angular 
frequency u is incident upon this filament, it will oscillate with a velocity that is 
proportional to the local flow velocity^]: 

Xi = C ij v° j e i & - jt - vot) . (3) 

and the detrmination of the coefficients will be the subject of subsequent 
computation. 

Whenever many vortex filaments are present and their mean separation is 
small by comparison with the sound wavelength they can be described in terms 
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of a continuous vorticity distribution u)(y, t) through the substitutions 

TjdaXr — Jtfyujfat) 

x — ► y . (4) 

Using this, and substituting (|3]) into (0) leads to the following approximate ex- 
pression for the amplitude ratio A of wave scattered into frequency v to incident 
wave: 

v r (f, v) ie iu ^l C2 vii 2 

A = -^—^ — ~ 2m e jikrikomC jm oJk(q, Az/), (5) 

v s c 2 \r\ 

where 6 is the scattering angle, a tilde denotes Fourier transform in space and 
time, q = [vrjc^ — k ) is the momentum transfer, Au = v — u is the frequency 
shift and 

2 ' jm 2 

Now, there are three important points that must be dealt with in order to turn 
the approximate expression (|^) into an exact one: The first is that the oscillatory 
vortex motion will generate not only second, but also first sound. The second 
is that second sound does not obey Eqn. (|l|), and the third is that, in addition 
to the local interaction between sound and vorticity just considered, there will 
be a nonlocal interaction due to the nonlinear coupling between the sound and 



Cjm (1 „ )3jm „ €jpm^p- 



vortical motions |15[], as well as additional effects due to mutual friction. 

In order to deal with the first point note that first sound will be radiated at 
wavelengths ~ c\jv and that, as a consequence the ratio of energy radiated into 
first sound waves to that radiated into second sound waves is of order c^/cf. We 
now choose to consider a situation near the lambda point, where this ratio is very 
small and we can neglect the first-to-second sound scattering. 

The second point is that second sound obeys the equation 

^-c^=0, (6) 

where w = v n — v s is the difference between the normal and superfluid velocities. 
However, in a scattering problem we wish to compare asymptotic states that are 
plane waves (globally for the incident wave, locally for the scattered wave), for 
which the relation 

PsO - 

V n = V s 

PnO 

holds. Here p s0 (resp. p n0 ) is the undisturbed superfluid (resp. normal) density. 
Since we have already agreed to be very near the lambda point, 

V n < V s (7) 
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and the superfluid velocity obeys the wave equation ([[]). Consequently, the am- 
plitude A given by (H) really gives the amplitude of second sound scattered by 
quantum vorticity, at least as far as the local interaction is concerned. 

Since temperature fluctuations are linearly related to superfluid velocity for a 
second sound wave, A will also give the ratio of incident to scattered temperature 
amplitudes. For the situation we have in mind, C2 ~ 10 3 cm/sec. One requirement 
that must be met in order for (|5]) to be valid is that the frequency of the incident 
wave must be high compared to the time scales of vortex flow on the length 
scale of the wavelength of the sound. Working in the 1-100 KHz range gives 
wavelengths in the 1 — 10~ 2 cm range, for which superflow velocities evolve in the 
10 _3 -10 _1 s range, so this requirement is amply met. 

The computation of the scattering of first sound is entirely analogous, it suf- 
fices to replace c<i by c±, the speed of first sound, in (|5]). This means that near 
the lambda point the scattered second sound can easily be three to four orders 
of magnitude larger than the scattered first sound. 

Dealing with the third point noted above, nonlocal interaction between sound 
and vorticity, as well as mutual friction, is harder. The proper way to do this is 



to start with the full nonlinear two- fluid equations |L4 



dV x 



dt 



+ V s • VK 



djpS) 
dt 



+ V ■ (pSV n 






-V/X + / 




(9) 
(10) 

fin 



with 



J 



n 



dp 
f 



PsV s + p n V n 

pSij + PnV n iV n j + p s V si V sj 

\ dp -SdT-^dW 2 
P 2p 

bpn -Q s A {Cj s A W) 



2p 



B'p r 

2p 



■oj q A W 



where p, p, S and T are mass density, pressure, entropy per unit mass and tem- 
perature respectively, lj s = V A V s and W = V n — V s . Thermal as well as viscous 
diffusivity effects have been neglected assuming the sound time scale to be small 
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enough. Terms describing the vortex self-induction have been omitted because 
they do not couple to the sound wave. Cubic terms in the entropy equation have 
been omitted because they will not contribute in the Born approximation to be 
considered shortly. 

The next step is to describe the flow as the sum of a fast "sound" (v) and 
slow "vortical" (u) parts: 

V = V + U !)<CM. 

That is, v describes the flow consisting of small deviations of density, temper- 
ature, entropy and pressure away from equilibrium values, around which the 
equations are linearized. These small deviations include both first and second 
sound. The vortical component u is a solution of the equations obtained by con- 
sidering constant density, with vortex filaments moving under their mutual and 
self-induction, possible external flows and appropriate boundary conditions. The 
interaction between the two modes is obtained by linearizing the equations not 
around a static value but around the vortical solution u. The result is 

|^-VV = F X (12) 



PnO d 2 S> 

PsoS 2 dt 2 



-V 2 T' = F 2 (13) 



where Fi and F 2 are interaction terms, quadratic in the dynamical variables, 
but linear both in u and in the sound variables that are described by primed 
quantities, with equilibrium values carrying an index "0" . This set of equations 
has the structure of a wave equation with source terms. If the "sources" F\ and 
F 2 are discarded, the usual equations describing sound (both first and second) 
are obtained. The solution to Eqn. (12) can be written as 



SOUND = SOUND mc + G*F (14) 

where SOUND is a two-vector describing first and second sound, G is a two-by- 
two matrix, the Green's function for the left-hand-side of (|T2"D, that can be found 
with standard Fourier techniques, and F is a two-vector with components Fi and 
F 2 . SOUND mc is a solution of the homogeneous equation. In an experiment such 
as we have in mind, SOU N D- inc is a plane, pure second sound wave, the incident 
wave, and the asymptotic behaviour of SOUND for large distances gives the 
scattered wave. It has both first and second sound components. Near the lambda 
point, c 2 <C c\ and p s0 <C p n0 . With these ingredients (|i~4"l) can be solved in a first 
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Born approximation that gives for the scattered temperature in the frequency 
domain T'(x, z/) scatt , 



TLtt ie lv ^ lc ^ 2 ( f . . B 



where 



n c\\r\ 

H(f)\ = n _ — \ 1 TqS dp' s 

K ' ~ K 2 } 2(cos^-i)Cp n0 ar'' 

C is the heat capacity, and 

*Hj = LOOij. 

Eqn. (|i5|) is our main result, for which a first approximation was provided in Eq. 
Taking into account the temperature dependence of the superfluid density p s 



and of the heat capacity near J\||16||, the last term in the right-hand side of H{6) 
can be estimated to be of order 2-3 at r = (T A — T)/T\ w 10 -3 , and to grow as 
r o.37 fj^jg value of r is not yet close enough to the A point for Eqns. (|^) to loose 
their validity. 

Let us now work out a few simple examples. First, one stationary, straight, 
vortex filament will have vorticity given by 

u(y,t)=T5^ 2 \y)y 3 

whose Fourier transform is, for 93 = 0, which is the only relevant situation in 
which the scattering plane is perpendicular to the filament, 

TL 

where L is the length of the filament or, in the case of a very long filament, the 
portion that is illuminated by the incident sound. Substitution into (|5|) shows 
that there is no frequency shift, as was to be expected for a static target. Using 
(H) and the frequencies we have quoted gives a typical scattered amplitude of 
10~ 5 — 10~ 7 , the higher amplitude corresponding to the higher frequency. This 
means that with a superconducting bolometer having a temperature resolution of 
10 nK it will be possible to detect the scattered second sound wave from a single 
quantum vortex if the amplitude of the incident wave is higher than 1 mK. These 
are certainly experimentally realizable conditions. 
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Consider next a circular vortex ring moving with velocity V. In general, the 
vorticity carried by a filamentary vortex X(a,t) is given by 

r r)X 

and, in the case of a circle, a convenient parametrization is given by 

X(a,t) = (R cos a, R sin a, Vt). 

Take the scattering plane to be the y\ — y 3 plane. In this case 

v 

qi = — cos 6 — k 

c 2 

q 2 = 

q 3 = — sin# (16) 
c 2 



so that a>2, the only relevant component, is given by 

PR V /•27T 

£o 2 = ——5(Av usin9) / da cos ffe -<*o*(«**-i)«» ^ . 

(27r) d c 2 Jo 

A Doppler shift in frequency arises from the vortex ring steady motion at velocity 
V. For rings whose radii are much smaller that the wavelength, koR <C 1 and 

—iTk R 2 . n <N . /A V . 
~ ^ — (cos 6* — l)()(Az/ z/sm#) 

87T^ c 2 

so that 

A ~ lO-^fcoi?) 2 , 

which appears to be very small even at 100 KHz and for rings of radius R ~ 10~ 4 
cm. Suppose however that we consider iV parallel vortex rings separated by a 
distance d, as might be the case for rings nucleating at a moving obstacle. In 
this case the scattering amplitude is enhanced by a factor iV along directions 
satisfying the Bragg condition 

dX sin 9 = n n = 1, 2, . . . 

Finally, attenuation techniques have suggested that the quantum turbulent 
state [17] may be characterized by an anisotropic vortex tangle with a difference 
in vorticity parallel and perpendicular to the direction of counterflow given ap- 
proximately by |0] 

Auj = TaA 
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where a is the anisotropy coefficient and A the vortex line density. Since the 
anisotropy is of order 1 and the experimental vortex line density can be as high 
as 10000, the corresponding scattering amplitude can be expected to be up to 
four orders of magnitude higher than that of the single quantum vortex. Also, 
in this case it is to be expected that scattering by the normal fluid vorticity will 
also be an important effect. 

To conclude, we have calculated the scattering of second sound by quantized 
vortices in superfluid Helium. Near the lambda point, this is an effect several 
orders of magnitude larger than first sound scattering and suggests a nonintrusive 
way of probing quantum vorticity. 
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